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We report on the fabrication aspects and calibration of the first large active mass (~ 15 g) mod- 
ules of SIMPLE, a search for particle dark matter using Superheated Droplet Detectors (SDDs). 
While still limited by the statistical uncertainty of the small data sample on hand, the first weeks 
of operation in the new underground laboratory of Rustrel-Pays d'Apt already provide a sensi- 
tivity to axially-coupled Weakly Interacting Massive Particles (WIMPs) competitive with leading 
experiments, confirming SDDs as a convenient, low-cost alternative for WIMP detection. 
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The rupture of metastability by radiation has been his- 
torically exploited as a method for particle detection. 
Perhaps its most successful application is the Bubble 
Chamber, where ionizing particles deposit enough local 
energy in a superheated liquid to produce vaporization 
along their wake. Apfel extended this concept in the form 
of Superheated Droplet Detectors (SDDs, a.k.a. Bub- 
ble Detectors), in which small drops (radius ~ 10 /xm) of 
the liquid are uniformly dispersed in a gel or viscoelas- 
tic medium. In a SDD the gel matrix isolates the fragile 
metastable system from vibrations and convection cur- 
rents, while the smooth liquid-liquid interfaces impede 
the continuous triggering on surface impurities that oc- 
curs in bubble chambers. The lifetime of the superheated 
state is extended, allowing for new applications: SDDs 
are increasingly popular as neutron dosimeters, where 
the nucleated visible bubbles provide a reading of the ra- 
diation exposure. SIMPLE (Superheated Instrument for 
Massive ParticLE searches) aims to detect particle dark 
matter using SDDs. We report here on the sensitivity at- 
tained at the early prototype stage, already comparable 
to the best achieved with competing technologies. 

In the moderately superheated industrial refrigerants 
used in SDDs, bubbles are produced only by particles 
having elevated stopping powers {dE/dx ^ 200 keV/nm) 
such as nuclear recoils. This is understood in the frame- 
work of the "thermal spike" model ||] , common to bubble 
chambers: for the transition to occur, a vapor nucleus of 
radius > Tc must be created, while only the energy de- 
posited along a distance comparable to this critical radius 
Tc is available for its formation. Hence, a double thresh- 
old is imposed: the deposited energy E must be larger 
than the work of formation of the critical nucleus, Ec, 
and this energy must be lost over a distance 0(rc), i.e., 
a minimum dE/dx is required. More formally P,H|: 



E> Ec = 47rr27/3e 
dE/dx > Ec/aVc, 



where rc — l 7/AP, 7(T) is the surface tension, AP = 
Pv—P, Pv{T) is the vapor pressure, P and T are the 
operating pressure and temperature, e varies in the range 
[0.02,0.06] for different liquids |,|, and a(T)~0(l) [|. 

Both thresholds can be tuned by changing the oper- 
ating conditions: keV nuclear recoils like those expected 
from scattering of WIMPs (currently the favored galactic 
dark matter candidates 0) are detectable at room T and 
atmospheric P, allowing for a low-cost search free of the 
complications associated to cryogenic equipment. Most 
importantly, the threshold in dE/dx provides an insen- 
sitivity to minimum-ionizing backgrounds that hamper 
the numerous WIMP detection efforts ||^. A mere <10 
WIMP recoils/kg target/day are expected and hence the 
importance of background reduction and/or rejection. 
SDDs of active mass 0(l)kg can in principle consider- 
ably extend the present experimental sensitivity 

Prompted by the modest active mass of commercially 
available SDDs (~0.03 g refrigerant/dosimeter) and the 
need to control the fabrication process, we developed 
a 80 60 bar pressure reactor dedicated to large-mass 
SDD production. It houses a variable-speed magnetic 
stirrer, heating and cooling elements and micropumps 
for catalyst addition (we nevertheless favored thermally- 
reversible food gels due to safety concerns in the han- 
dling of synthetic monomers). The fabrication of 1 I SDD 
modules containing up to 3% in superheated liquid starts 
with the preparation of a suitable gel matrix; ingredients 
are selected and processed in order to avoid alpha emit- 
ters, the only internal radioemitters of concern A 
precise density matching between matrix and refrigerant 
is needed to obtain a uniform droplet dispersion, making 
water-based gels inadequate unless large fractions of inor- 
ganic salts are added, which can unbalance the chemistry 
of the composite and contribute an undesirable concen- 
tration of these contaminants . We find that glycerol 
is for this and other reasons an additive of choice. It is 
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purified using a bed of pre-eluted ion-exchanging resin 
specifically targeted at actinide removal. Polymer addi- 
tives and gelating agent are washed in a resin bath. All 
components are forced through 0.2 /im filters to remove 
motes that can act as nucleation centers. The resulting 
mixture is outgassed and maintained above its gelation 
temperature in the reactor. The refrigerant is distilled 
and incorporated to this solution at P >> Pv{T) to 
avoid boiling during the ensuing vigorous stirring. After 
a homogenized dispersion (r ~ 30 ± 15 /im) of droplets is 
obtained, cooling, setting and step- wise adiabatic decom- 
pression produce a delicate entanglement of superheated 
liquid and thermally-reversible gel, the SDD. The detec- 
tors are refrigerated and pressurized during storage to 
inhibit their response to environmental neutrons. 

SDDs can bypass the listed problems associated to 
a former bubble chamber WIMP search proposal, 
but are not devoid of their own idiosyncrasies. For in- 
stance, the solubility of hydrogen-free refrigerant liquids 
in water-based gels is small (e.g., 0.002 mol/kg bar for 
R-12, CCI2F2), yet sufficient to produce unchecked bub- 
ble growth via permeation after few days of continuous 
SDD operation. The engorged bubbles lead to fractures, 
spurious nucleations and depletion of the superheated 
liquid (commercial gel-based SDDs are designed for few 
hours of exposure before recompression [ p^ , a cycle that 
can be repeated a limited number of times). To achieve 
the long-term SDD stability needed for a WIMP search, 
we employ a multiple strategy: fracture formation can 
be delayed under a moderate P ^ 2 atm, or by choos- 
ing refrigerants with the lowest solubility in the matrix 
(~ 0.0003 mol/kg bar for R-115, CzClFr,). Structure- 
making inorganic salts produce a "salting-out" effect, 
i.e., further reduce the refrigerant solubility. Their use 
being inadvisable for the reasons above, we introduce 
instead polymers known to have a similar effect [p^ , 
such as polyvinylpyrrolidone (PVP). As a result of these 
measures, present SIMPLE modules are stable over ~40 
d of continuous exposure. Another example of SDD- 
specific problems is the formation of clathrate-hydrates 
on droplet boundaries during fabrication or recompres- 
sion. These metastable ice-like structures are inclusions 
of refrigerant molecules into water cages that shorten 
the lifetime of superheated drops encrusted by them via 
transfer mechanisms still not well understood |T5| . Their 
presence may be responsible for a long-lived spurious nu- 
cleation rate observed in R-12 SDDs following fabrication 
[0 . This is addressed in SIMPLE with the addition of 
polymers such as PVCap or PVP, which act as kinetic 
inhibitors in their growth JTif , and by use of large molec- 
ular size refrigerants like R-115, for which the formation 
of most hydrates is stoichiometrically forbidden |l^,|l^ . 

Prototype modules are tested in an underground 
gallery. The 27 m rock overburden and ^ 30 cm paraf- 
fin shielding reduce the flux of muon-induced and cos- 
mic fast neutrons, the main source of nucleations above 



ground. Inside the shielding, a water-|-glycol thermally- 
regulated bath maintains T constant to within 0.1°C. 
The characteristic violent sound pulse accompanying va- 
porization in superheated liquids Jl7| , p^ is picked-up by a 
small piezoelectric transducer in the interior of the mod- 
ule, amplified and stored. Special precautions are taken 
against acoustic and seismic noise. Fig. 1 displays the 
decrease in spontaneous bubble nucleation rate brought 
by progressive purification of the modules. 
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FIG. 1. SDD background at 90 m.w.e. and P=2 atm, fol- 
lowing cumulative steps of cleansing; histogram: double distil- 
lation of water and microfiltration, •: distillation of refrigerant 
and glycerin purification, o: gelatine and PVP purification. In- 
sert: Rate induced in a calibration (0.2 Bq ^''^Am/g gel). The 
theoretical values of Tar and Ta are indicated (the dashed line 
is an numerical simulation containing no free parameters). 

The response of smaller SDDs to various neutron 
fields has been extensively studied |^,|l^ and found 
to match theoretical expectations. However, large-size, 
opaque SDDs require independent calibration: acoustic 
detection of the explosion of the smallest or most distant 
droplets is not a priori guaranteed. The energy released 
as sound varies as (Pv - Pf'"^ [Hi, making these addi- 
tional characterizations even more imperative for SDDs 
operated under P> \ atm. Two separate types of calibra- 
tion have been performed to determine the target mass 
effectively monitored in SIMPLE modules and to check 
the calculation of the T, P-dependent threshold energy 
Ethr above which WIMP recoils can induce nucleations 
(defined as the lowest energy meeting both conditions in 
Eq. (1) [H^JH). First, a liquid Am source (an alpha 
emitter) is diluted into the matrix while still in the so- 
lution state. Following Eq. (1), the 5.5 MeV alphas and 
91 keV recoiling ^'^^Np daughters cannot induce nucle- 
ations at temperatures below Ta and Tar, respectively 
Q. The expression a = 4.3 (pu/p/)^^'^ where Pv{T), 
pi (T) are the vapor- and liquid-phase densities of the re- 
frigerant, correctly predicts the observed Tq, for both R- 
12 and R-115 at P=l and 2 atm. In the same conditions, 
the theoretical value of e |2|] for these liquids (e ^0.026, 
neglecting a small T, P dependence) generates a good 
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agreement with the experimental Tar (Fig. 1, insert). 
Prior to extensive component purification, the spectrum 
in non-calibration runs (Fig. 1, histogram) bears close 
resemblance to that produced by ^^^Am spiking (Fig. 1, 
insert); the initial presence of a small (~ 10"'* pCi/g) 
^^®Th contamination, compatible with the observed rate, 
was confirmed via low-level alpha spectroscopy. Three 
regimes of background dominance are therefore delimited 
by vertical lines in Fig. 1: the sudden rise at T ^ 15°C 
originates in high-dE / dx Auger electron cascades follow- 
ing interactions of environmental gammas with CI atoms 
in the refrigerant Q,^- The calculated Ec for R-115 at 
T = 15.5°C andP=2 atm is 2.9 keV, coincidental with the 
binding energy of K-shell electrons in CI, 2.8 keV (i.e., 
the maximum E deposited via this mechanism). Thus, 
the onset of gamma sensitivity provides a welcome addi- 
tional check of the threshold in the few keV region. 

Alpha calibrations are not suitable for a rigorous deter- 
mination of the overall sound detection efficiency because 
a large fraction of the added emitters drifts to gel-droplet 
boundaries during fabrication, an effect explained by the 
polarity of actinide complex ions ]2^ and dependent on 
matrix composition. While this migration does not affect 
Ta nor Tar, it enhances the overall nucleation efficiency 
in a somewhat unpredictable manner ]2^ . To make up for 
this deficiency, SIMPLE modules have been exposed to a 
^^^Cf neutron source at the TIS/RP calibration facility 
(CERN). The resulting spectrum of neutron-induced flu- 
orine recoils (Fig. 2, insert) mimics a typically expected 
one from WIMP interactions. A complete MCNP4a ||2|] 
simulation of the calibration setup takes into account the 
contribution from albedo and thermal neutrons. The ex- 
pected nucleation rate as a function of T is calculated 
as in 1^,^: cross sections for the elastic, inelastic, (n,Q!) 
and (n,p) channels of the refrigerant constituents are ex- 
tracted from ENDFB-VI libraries. Look-up tables of the 
distribution of deposited energies as a function of neutron 
energy are built from the SPECTER code ^ , stopping 
powers of the recoiling species are taken from SRIM98 
p7| . Since T was continuously ramped up during the ir- 
radiations at a relatively fast l.l°C/hr, a small correction 
to it (< 1°C) is numerically computed and applied to ac- 
count for the slow thermalization of the module. Depend- 
ing on r, the value of Ethr for elastic recoils in fluorine 
(the dominant nucleation mechanism in R-115) is set by 
either condition in Eq. (1), the other being always ful- 
filled for E > Ethr ||9|,p0||. The handover from the second 
to the first condition at T above - 5.5°C (^ 2.5°C) for 
P =2 atm {P —1 atm) is clearly observed in the data 
as two different regimes of nucleation rate (Fig. 2). A 
larger-than-expected response, already noticed in R-12 
1^, is evident at low T: the calculated Ethr there is too 
conservative (too high). This behavior appears well be- 
low the normal regime of SDD operation (which is at T 
high enough to have Ethr — Ec) and therefore does not 
interfere with neutron or WIMP detection. However, it 



is interesting in that it points at a higher than normal 
bubble nucleation efficiency from heavy particles, as dis- 
cussed in early bubble chamber work A best-fit to 
the overall normalization of the Monte Carlo over the full 
data set (Fig. 2, dotted fines) yields the fraction of refrig- 
erant mass monitored with the present sound acquisition 
chain, 34 ± 2% (74 ± 4%) of the total at P=2 atm (P=l 
atm), a decisive datum to obtain dark matter limits. 
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FIG. 2. ^^^Cf neutron calibration of SIMPLE modules at 
the TIS/RP bench (CERN), compared with Monte Carlo ex- 
pectations (dotted lines, see text). The signal-to-noise ratio 
was > 30 at all times. Insert: calculated energy spectrum of 
F recoils during the irradiations. 

The installation 500 m underground of modules iden- 
tical in preparation and sound detection system to those 
utilized in ^^^Cf calibrations started in July 1999. A de- 
commissioned nuclear missile control center has been con- 
verted into an underground laboratory facilitating 
this and other initiatives. The characteristics of this site 
(microphonic silence, unique electromagnetic shielding 
|28| ) make it specially adequate for rare-event searches. 
Modules are placed inside a thermally-regulated water 
bath, surrounded by three layers of sound and thermal 
insulation. A 700 1 water neutron moderator, resting on a 
vibration absorber, completes the shielding. Events in the 
modules and in external microphones are time-tagged, al- 
lowing to filter-out the small fraction (~ 15%) of signals 
correlated to human activity in the immediate vicinity of 
the experiment. P and T are continually logged. The sig- 
nal waveforms are digitally stored, but no event rejection 
based on pulse-shape considerations [ p^ is performed at 
this stage, eluding the criticisms associated to some 
WIMP searches in which large data cuts are made. 

The raw counting rate from the first SIMPLE module 
operated in these conditions appears in Fig. 3. Account- 
ing for sound detection efficiency and a 62% fluorine mass 
fraction in R-115, limits can be extracted on the spin- 
dependent WIMP-proton cross section gyi/p (Fig. 3). The 
cosmological parameters and method in [ pO[ are used in 
the calculation of WIMP elastic scattering rates, which 
are then compared to the observed uncut nucleation rate 
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at T = 10°C or 14°C, depending on WIMP mass. The 
expected nucleation rate at T (i.e., integrated for re- 
coil energies > Ethr{T)) from a candidate at the edge 
of the sensitivity of the leading DAMA experiment ||] 
(~ 1.5 • lO'* kg-day of Nal) is offered as a reference in 
Fig. 3: SIMPLE sensitivity is presently limited by the 
large statistical uncertainty associated to a short expo- 
sure, and not yet by background rate. A considerable 
improvement is expected after the ongoing expansion of 
the bath to accommodate up to 16 modules. In paral- 
lel to this, plastic module caps are being replaced by a 
sturdier design: runs using refrigerant-free modules show 
that a majority of the recorded events arise from pressure 
microleaks, correlated to the sense of T ramping, able to 
stimulate the piezoelectric sensor. It must also be kept in 
mind that a T— independent, flat background implies a 
null WIMP signal, albeit this eventual approach to data 
analysis can only be exploited after a large reduction in 
statistical uncertainty is achieved. 
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FIG. 3. 95% C.L. limits on awp extracted from 0.19 
kg-day of SDD exposure, compared with other experiments 
[8]. "MSSM" marks the tip of the region where a lightest su- 
persymmetric partner is expected. Insert: Counting rate in the 
module (9.2 ± 0.1 g R-115, AT = -0.75°C/day). The dotted 
line is the expected signal (corrected for 34% detection effi- 
ciency) from a WIMP of mass =10 GeV and awp = 5 pb: 
present sensitivity is still limited by low statistics. 

The importance of the spin-dependent WIMP inter- 
action channel (where F is the optimal target |^l|) has 
been recently stressed by its relative insensitivity to CP- 
violation parameter values, which may otherwise severely 
reduce coherent interaction rates Nevertheless, 
CF^Br modules able to exploit coherent couplings are 
presently under development. The intrinsic insensitivity 
of SDDs to most undesirable backgrounds, low cost of 
materials involved and simplicity of production and op- 
eration opens a new door to dark matter detection. 
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